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USER’S GUIDE TO BREAK], The Berkeley Algorithm for 
Breaking Window Glass in a Compartment Fire 
Aruna A. Joshi and Patrick J. Pagni 
Mechanical Engineering Department 


University of California at Berkeley 
Berkeley, CA 94720 


Abstract 


This is an instructional manual for the computer program BREAKI], the 
Berkeley Alogorithm for Breaking Window Glass in a Compartment Fire, version 1.0. The 
thermal response of window glass exposed to a compartment fire and its time to breakage are 
calculated. The following glass properties, treated here as constant, are required as input: 
thermal conductivity, thermal diffusivity, absorption length, breaking stress, Young’s modulus, 
coefficient of linear thermal expansion and emissivity. In addition, the ambient temperature, 
the hot layer temperature, heat transfer coefficients and emissivities and the direct flame radia- 
tion histories are needed as input. The window geometry is also required. The output 
describes the evolution of the window temperature field, T(x,t), where x is the direction nor- 
mal to the glass pane, culminating in the window breakage time, t,. This manual contains a 
summary of the mathematical model, sample thermophysical and mechanical data tables, some 


worked examples and advice on using the program. 


i) 


1. INTRODUCTION 


Window glass breaking plays a key role in fire growth. Before breaking, the window 
acts as a wall and after breaking, as a vent. It is therefore necessary to know the time at 
which a window breaks in order to model a compartment fire. Windows break in a fire due to 
tensile stresses created in the portion of the glass shaded by the frame as the exposed central 
region heats up and expands [1]. Since glass has a low tensile strength, a reasonably small 
stress ( ~ 40 MPa) is enough to cause breakage. In FIRST[2], the code for fire simulation in 
a room, there is a provision for opening a vent (i.e. a window) at a specified time. This pro- 
gram provides the specified time, i.e., when the window breaks and becomes a vent, for any 
user specified flame radiation and layer temperature histories. The units for all the input and 


Output are SI. 


This document provides descriptions of the program BREAKI. The next section 
discusses the mathematical model and computer program. The third section describes the 


required input. The fourth section gives example sessions and the fifth explains the output. 


2. MATHEMATICAL MODEL 


The computer program BREAK I provides a solution to the transient problem of heating 
a window by exposing it to a fire and its hot products. The window glass absorbs direct radia- 
tion from the flame within its depth and is heated by radiation and convection from the sur- 
rounding hot gases while it loses heat to the cooler outside ambient atmosphere. In the shaded 
portion of the window glass within the frame, large tensile stresses develop. These tensile 
stresses are directly proportional to the average temperature difference between the shaded and 
the unshaded portion of the window glass. Figure 1 shows a schematic diagram of window 
glass exposed to a fire. L is the glass thickness, s is the evils of the window shaded by the 


frame and H is the half-width or half-length of the window. 


In a review of fire physics [3] at the Second International Symposium on Fire Safety Sci- 


ence, a simple strain criterion for the glass temperature increase AT required to break windows 


rm > 


Sec. A-A 


Fig. 1. Window Geometry 


Sf, force balance factor 
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Fig. 2. Increase in temperature rise at breaking due to compres- 
sion in the heated glass. For a step temperature profile, 
f = 1+s/H, for a hyperbolic tangent temperature profile, 
f = 2/[tanh(s/L) + L/(s+H )(In((cosh(H/L )icosh(s/L )))). 


in fires, 
BAT =O, 72, (1) 


was suggested, where B is the thermal coefficient of linear expansion, 6, is the tensile strength 


at breakage and E is the Young’s Modulus of the glass. 


This relation was recently improved [4] to include a multiplying factor g(~1) on the nght 
hand side of Eq. (1). to account for compression in the unshaded glass. The values for g in 
Fig. 2 are obtained from an over-all force balance in the plane of the pane. For large shading, 
s/L22, and fast heating, at/s*<1, where a = k/pc is the thermal diffusivity of glass, the 
shaded edge of the window glass remains quite near the initial temperature, T;. Hence, the 
temperature change AT at which breakage occurs is the increase in the average temperature of 
the central glass portion given by 


if 
AF = (UL)|T (x t)dx - T; = g5,/EB. (2) 
0 


Thus the goal here is to calculate the temperature in the large central section of the glass 
as a function of depth into the glass, x, and time, t. The time to initial window breakage is 
determined. The question of additional time required for all of the glass pieces to fall from the 


window frame [5] is not addressed. The governing equation for this system is 


2 | —x/l 
sya an cena AO ret (3) 


ot ox : 


where /(t) is the incident radiative flux directly from the fire which is at sufficiently short 
wavelengths that its distributed internal absorption needs to be included [6] and / is the decay 


length in the glass [7]. Radiation from the hot layer is included in the following boundary 


conditions 
t= 0; f= J ge. (4) 
bli tesa 56 4 4 
Pee) iy ret = hy(t)(T s(t )—-T (0,t)) + €2..(¢ OT 3. (t) -eoT (0,1) = g(t), (5) 


te Le oT =p (T(L ag )\-T 20) uF eoT (L Ly +s ee er = i(¢). (6) 
Bf 


Here, & is the thermal conductivity, p is the density, c is the specific heat, h(t) and A, are the 
heat transfer coefficients on the exposed and unexposed sides respectively. T>..(¢) is the com- 
partment hot layer temperature and 7. is the ambient temperature, o is the Stefan-Boltzmann 


constant, €, €;,, and €5,, are the emissivities of glass, ambient and hot layer respectively. 


The equations are nondimensionalized using the following variables 


(7) 


The equations are solved using Laplace Transforms and Newton’s method. An advantage of 
this technique [8], suggested by H. R. Baum, is that the inversions need be performed only at 


the boundaries for the surface temperatures. 


We showed in Ref. 9 (Fig. 3 (steady parameter solution)) that the linear interpolation of 
temperature to obtain the average temperature underestimates the time to breakage. However, 
the linearity increases as T increases. Then, a reasonable approximation for average tempera- 
ture would be to obtain a function 


l 
Jog.naé 


X(t) = (8) 


0 
@(0,t) + O(1,T) 


and approximate 6(E,t) = ~(0(0,t) + O(1,t)). In the limit of a linear temperature profile, 


x = 0.5. Figure 3 shows the variation of x with t for the case of Fig. 3 in Ref. 9 and the fit 


X= ~ orf (Vt) which appears to be reasonable. For tS0.05, the semi-infinite unsteady tem- 


perature profile applies, while for t>1, x asymptotically approaches 0.5. In reality, the value 
of x will be even smaller than that estimated from Eq. (8) i.e. the temperature profile will be 
even more nonlinear since unsteadiness usually increases heating causing a higher 9(0,7). 


Nevertheless, this variation of x would still be better than the linear value of 0.5. 
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BREAK] is written in FORTRAN 77. It has about 2200 lines of code including the 
comment statements. No description of the numerical procedure is provided (please refer to 
Pagni and Joshi [4]). However, an iterative procedure is utilized which checks the solutions by 
doubling or halving the time steps. If the difference is greater than the specified accuracy 
(default being (soln (Ar )—soln ((At /2)\/soln (Ar )S107*), the time step is halved again and the 
computations are repeated. Due to memory constraints, the number of time step reductions is 
restricted to 3 after which a message is displayed saying that the maximum allowable error is 


too small and the computations are stopped. 
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Figure 3. Variation of the multiplying factor y = TENE FU as a function of t for a 


constant parameter example of Fig. 3 of Ref. 9. 


3. INPUT 


The physical and mechanical properties and. geometry of glass, fire parameters and 
numerical parameters required by the program are listed in Table | along with the input format 


and a suggested range of properties [10-13]. 


The hot layer temperature, the flame radiation, the heat transfer coefficient on hot layer 
side and emissivity of hot layer functions are entered at user specified time increments. The 
program then interpolates between the end points and builds a function. After the last entered 


time, the program assumes that the value remains constant. 


There are two ways to input parameters. The first way is by changing default values and 
marking the altered file as an input file. The second one is by reading in data from a pre- 
existing file. Changes can be made to the input read from an input file using the following 


menu which enables the user to pick any desired area for alteration. 


1) Change Properties of Glass 

2) Change Geometry of Window 

3) Change Coefficients 

4) Change Flame Radiation 

5) Change Gas Temperature 

6) Change Heat Transfer Coefficient on Hot Layer Side 
7) Change Emissivity of Hot Layer 


8) Change Numerical Parameters 


The program has been designed so that in the interactive mode, it reads the data in free 
format. After each entry the "carriage return" or "enter" key has to be hit. After changes are 
made in a particular section, the program shows the changed data on the screen and gives the 


user another chance to make changes in that section. 


hi " 


When typing input, if an error is made, e.g. instead of "y" a "t" is typed, or instead of a 
character, a number or vice versa, is typed, the question is re-asked. Also, if a number outside 


the permitted range is typed, the question is re-asked. 


Table 1. Input 


| | } 


| Name 
| 


| Thermal conductivity k 


; W/mK El2.4 | 0.6-1.3 | 


Thermal diffusivity | a m/s E12.4 | 3.0-4.0x1077 | 
| Absorption length | J m B12.4 | 0.9-1.5 xl03 | 
| Breaking stress ls aO;, N/m- E12.4 | 2.0-13.8 x10’ | 
Young's modulus | E Nim* | E12.4 -7.0-7.25 x10 
| Thermal coefficient of linear expansion | B oe! | E12.4 | 8.5-9.5 xl0~° | 
| Thickness of the glass | L m | F8.4 | 0:123.0 x1 | 
| Shading thickness ls m F840 | 0.3-5 x10? 
"-Half-width of the window |A m | F8.4 | 0.05-5.0 | 
| Heat transfer coefficient on exposed side hx{t) Wim?K | F7.2 | 5.-80. | 


| Heat transfer coefficient on unexposed side | h, 


| Ambient temperature on unexposed side 


| | 


270.-320. 


| Emissivity of glass ae 


\ ' | | | 
SE a a ee 


_ Emissivity of ambient 


| 
{ 
| 
i 


| Emissivity of hot layer 


F6.2 


i \ 1 
eee eee a 


0. Soin 
| Onan 


Odett 


"Time tf s F12.2_ | 0.0-1500. 
"Flame radiation flux | I(t) W/m? | F12.2. | 0.0-10000. 
Hot layer temperature | Toate K | F12.2. | 300.-1300. | 
| 
s Oa oe ee be | 


| Time increment At 


| Maximum error allowed in solutions 


| Maximum run time | t 


| Output time interval 


| &e 


F820 | 
cana = ae Tee ee 


| F72 


100.-1500. 


| 


| 1-100. | 


4. EXAMPLE SESSIONS IN INTERACTIVE MODE 


This section shows several example sessions using different options. Execute a run by 
typing BREAKI1. All the responses are either y, standing for yes, n, standing for no, numbers 
or a filename. End your entry by typing any two letters. All questions must be answered with 
some response or the computer will continue to wait for a response. Comments in parentheses 
do not appear on the screen; they are included to explain the options. Responses by the user 
are underlined. The input for the default example 4.1, was obtained by modeling a particular 
compartment fire with the program FIRST [9]. The data are shown in the output statement 
given in section 5. The compartment was large and the fire was far from the window so that 
I(t) =0. The fire grew exponentially with a 30s doubling time. In example 4.2, you will 
interact with the nn by changing a single input parameter. Example 4.3 is similar, but 
here you change an entire input field. In example 4.4, you recall the input file saved in exam- 
ple 4.1, modify it and use it as input. You have three options for changing values in input 
files: 1. using such files unchanged, 2. editing the input file with a separate editor and 3. edit- 


ing the input file interactively. The latter technique is shown in example 4.4. 


4.1. Default Run 


This is the shortest session which is the default case. You will save your input as a file 


named exl.in to be modified in example 4.4. 


BREAKS VERSION 1.0 
A- A. JOShi and P. J. Pagni 
Department of Mechanical Engineering 
University of California at Berkeley 
Berkeley, CA 94720 


This program calculates the temperature history of the surfaces 
of a glass window for given fire parameters. 


The calculations are stopped when the glass breaks. 


Bit RETURN) to continue 


KKK HHH KKK KKK KKK KK KKK KKK KKK KKK KKK KKK RK KKK KKK KK KKK KKRKEKKK 
* Warning: Read the disclaimer before using this program 

~~ awww wm wm em wm em em wm e@ ew ewe wm @ @ @ MP @@ @@ ewe @ ZT @@ @ 2 eee @ MO ew @ @O 2m ee ee eee oe — 

* This program was developed under grants from the 

* National Institute of Standards and Technology, however 
* no certification or applicability is implied. Many 

* assumptions are made concerning the physical properties 
* of glass and the phenomenology of breakage. Only you are* 
* responsible for your calculations with this program. Use* 
* should be restricted to predicting when ordinary glass * 
* in ordinary windows may break under the assumptions in * 
* your choice of input parameters. Do not use the 
* 

*x 

« 

* 

* 

* 

* 


+ 4+ + + + 


program to calculate conditions under which glass will 
not break. The properties of a particular piece of 
glass and the accuracy of each set of input are both 
sufficiently uncertain that no confidence can be placed 
in any calculation which predicts that a window will 


not break. 
KEKE KEKE KKK KKKKKE KEKE KKEKEKRRKEKKRKEKEKRKEKKREKKEKKKKRKKKK 


* 22% AD 1. 2 


Pleaseiiwrites(y’ for yes and ’n° forsno 
for every question. 


Do you want to use the default values for properties, geometry, 
fire and numerical parameters? (y/n) 
y (y calculates the temperature histories for default values) 
(n see example 4.2) 


Would you like to save your input data in a file?(y/n) 
vg 


What would you like to call it?(must be in DOS format) 
exl.in (your choice, I’ iil choose ¢2:12-1n) 


Would you like to save the output in a file? (y/n) 


What would you like to call it?(must be in DOS format) 
exl.out {your choice), Tollischoose.e=1 our) 

Begin run 

(See the example output in section 5) 
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1] 


4.2. Change Input Parameters 


In this run only the thickness of the glass is changed. The procedure for changing data 
in the areas COEFFICIENTS and NUMERICAL PARAMETERS is identical to the changes 


made in the area of GEOMETRY as in this run. 


Do you want to use the default values for properties, geometry, 

fire and numerical parameters? (y/n) 
n 

Do you want to change the default values interactively or read in from 
a pre-existing file? (y/n) 


NOTE: To read in values from a file, you have to use a file saved ina 
previous run and if needed change the data according to the format specified in 
Table I of the manual. 
y=change default values interactively, n=read from a file 
y (y calls the default values in the program for any changes to be made) 
(n reads data from a user specified file) 
To change the input in any of the following sections,enter the item number 
of your choice 


1)NO MORE CHANGES 

2) CHANGE PROPERTIES OF GLASS 

3) CHANGE GEOMETRY OF THE WINDOW 

4) CHANGE COEFFICIENTS 

5) CHANGE FLAME RADIATION 

6) CHANGE GAS TEMPERATURE 

7) CHANGE HEAT TRANSFER COEFF. ON HOT LAYER SIDE 
8) CHANGE EMISSIVITY, OF HOT, LAYER 

9) CHANGE NUMERICAL PARAMETERS 


3 


GEOMETRY 

1.Glass thickness {m]= 0.0063 
2.Shading thickness [m]J= 0.0150 
3.Half-width [m)= 0.5000 


Write item numbers and values you wish to change separated by comma. 

Hit the RETURN or ENTER key after entering each pair. End your entry series 
by typing any 2 characters. (e.g. 1,.00635 RETURN 2,.02 RETURN EN) 

1, 6 O05 
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GEOMETRY 

1.Glass thickness [m]= 0.0050 
2.Shading thickness [m])= 0.0150 
3.Half-width [m]= 0.5000 


Would you like to make any more changes in this section? (y/n) 
n (y asks to enter numbers again) 
= (n returns to the menu) 


To change the input in any of the following sections,enter the item number 
of your choice 


1)NO MORE CHANGES 

2) CHANGE PROPERTIES OF GLASS 

3) CHANGE GEOMETRY OF THE WINDOW 

4) CHANGE COEFFICIENTS 

5) CHANGE FLAME RADIATION 

6) CHANGE GAS TEMPERATURE 

7) CHANGE HEAT TRANSFER COEFF. ON HOT LAYER SIDE 
8) CHANGE EMISSIVITY OF HOT LAYER 

9) CHANGE NUMERICAL PARAMETERS 


Would you like to save your input data in a file? (y/n) 
ee 

What would you like to call it?(must be in DOS format) 
ex2.in (your choice, I’1l1 choose ex2-in) 


Would you like to save the output in a file? (y/n) 


What would you like to call it?(must be in DOS format) 
ex2.out (your choice, I° lI choose exzc0uc, 


Begin run 


—" 
ty) 


4.3. Change Input Fields 


This session describes the case of a fire next to the window. Hence, changes are made in 
the area FLAME RADIATION. In the default run the flux from the fire is equal to zero. If 
the run time is greater than the last entry, the flux stays constant until the maximum run ume. 
The procedure for changing the GAS TEMPERATURE, HEAT TRANSFER COEFFICIENT 
ON HOT LAYER SIDE and EMISSIVITY OF HOT LAYER is identical to that for changing 


FLAME RADIATION. 


To change the input in any of the following sections,enter the item number 
of your choice 

1)NO MORE CHANGES 

2) CHANGE PROPERTIES OF GLASS 

3) CHANGE GEOMETRY OF THE WINDOW 

4) CHANGE COEFFICIENTS 

5) CHANGE FLAME RADIATION 

6) CHANGE GAS TEMPERATURE 

7) CHANGE HEAT TRANSFER COEFF. ON HOT LAYER SIDE 
8) CHANGE EMISSIVITY OF HOT LAYER 

9) CHANGE NUMERICAL PARAMETERS 


FLAME RADIATION 


This program can handle time varying flame radiation flux. The flux 
is specified with values at discrete intervals of time. The program 
then linearly interpolates the flux between time intervals. 


Values for flux are 


point # time [s] Luce TWO. 2, 
. a (ep, 
ws 1000.00 On 


What type of change would you like to make? 
1. Add time/s 

Cothewer ale 

3. Enter new data 
3 
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Enter time in secs and flux in W/m*2 separated by a comma. 

Hit the "RETURN" or "ENTER" key after entering each pair. End your 
entry series by typing any 2 characters. (e.g. 0,10000. RETURN EN) 
0,0 


200,10000 
500,3000 
en 
point # time [s] flux (W/m*2] 
al OF OR 
7 200.00 1000.00 
3 500.00 3000.00 


Would you like to make any more changes in Spe section? (y/n) 
n (y asks to change values again) 
- (n returns to the menu) 


To change the input in any of the following sections,enter the item number 
of your choice 


1)NO MORE CHANGES 

2) CHANGE PROPERTIES OF GLASS 

3) CHANGE GEOMETRY OF THE WINDOW 

4) CHANGE COEFFICIENTS 

5) CHANGE FLAME RADIATION 

6) CHANGE GAS TEMPERATURE 

7) CHANGE HEAT TRANSFER COEFF. ON HOT LAYER SIDE 
8) CHANGE EMISSIVITY OF HOT LAYER 

9) CHANGE NUMERICAL PARAMETERS 


Would you like to save your input data in a file? (y/n) 
yy, 

What would you like to call it?(must be in DOS format) 
ex3.an (your choice, Ill enhoese ex3, 15) 

Would you like to save the output in a file?(y/n) 

ee 

What would you like to call it?(must be in DOS format) 
e=3,0u0t (your choice, I’ ll choose ex3 out) 


Begin run 


4.4. Input Files 


This sample uses the file saved in the previous run of example 4.1 as the input. The 
input file can either be edited by using the editor on the computer or by changing values 
interactively. In this case, numerical parameters and the gas temperature history are changed 


interactively. 


Do you want to use the default values for properties, 
fire and numerical parameters? (y/n) 
n 


geometry, 


Do you want to change the default values interactively or read in from 
a pre-existing file? (y/n) 


NOTE: To read in values from a file, you have to use a file saved ina 
previous run and if needed change the data according to the format specified in 
Table I of the manual. 


y=change default values interactively, n=read from a file 


Please write the file name 
exl.in 


(must be in DOS format) 


Would you like to change any input data? (y/n) 


To change the input in any of the following sections,enter the item number 
of your choice 


1)NO MORE CHANGES 


2) CHANGE PROPERTIES OF GLASS 

3) CHANGE GEOMETRY OF THE WINDOW 

4) CHANGE COEFFICIENTS 

5) CHANGE FLAME RADIATION 

6) CHANGE GAS TEMPERATURE 

7) CHANGE HEAT TRANSFER COEFF. ON HOT LAYER SIDE 
8) CHANGE EMISSiIVGLY Ob HOs yLAYER 

9) CHANGE NUMERICAL PARAMETERS 
9 

NUMERICAL PARAMETERS 

1.Maximum fractional error in soln= 0.000100 
2.Size of time step [s]= 1.000 

3.Maximum run time [s]= 250.00 

apcome interval. for output fst= 10.00 
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Write item numbers and values you wish to change separated by comma. 


Hit the RETURN or ENTER key after entering each pair. 


by typing any 2 characters (€.'9. 


1,.001 
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1,1.e-4 RETURN EN) 


1.Maximum fractional error in soln= 0.001000 


2.Size of time step [s]= 1.000 
3.Maximum run time [s]= 250.00 
4.Time interval for output [s]= 


Would you like to make any more 


£03.06 


End your entry series 


changes in this section? (y/n) 


To change the input in any of the following sections,enter the item number 


of your choice 


1)NO MORE CHANGES 

2) CHANGE PROPERTIES OF GLASS 

3) CHANGE GEOMETRY OF THE WINDOW 
4) CHANGE COEFFICIENTS 

5) CHANGE FLAME RADIATION 

6) CHANGE GAS TEMPERATURE 

7) CHANGE HEAT TRANSFER COEFF. ON 
8) CHANGE EMISSIVITY OF HOT LAYER 
9) CHANGE NUMERICAL PARAMETERS 


GAS TEMPERATURE 


HOT LAYER SIDE 


This program can handle time varying hot ambient temperature. The 
temperature is specified with values at discrete intervals of time. 
The program then linearly interpolates the temperature between time 


intervals. 
Values for hot layer temp are 
point # time [s! 
On 
10500 
20.00 
30.00 
40.00 
50.00 
60.00 
70.00 
80.00 
90'..00 
100.00 
110.00 
£20500 
140.00 
160.00 
180.00 
200.00 


PRP eP 
OWDIDUMPWNHKH 


NYHA UW PWD bE 


PPRPRPH 


temperature 
300.00 
303730 
303.. 96 
304.60 
SUSso7 
306.85 
308.20 
Bi0.09 
312.64 
Bhi Tele gS 
319356 
Ba tee dah she 
Salts 
353.234 
206,20 
437.97 
S31.07 


(K] 


What type of change would you like to make? 


1. Add time/s 

2. Alter data 

3. Enter new data 
ic) 


Enter time in [s] 


and temperature in [K] 


separated by a comma. 


Hit the “RETURN” or the "ENTER" key after entering each pair. 


End your entry series by typing any 2 characters. 


Pou 

10,400 

20,600 

30,700 

4054750 

en 

point # time {s] 

z QO. 
Z LTOe 00 
3 20800 
4 30.00 
5 40.00 


Would you like to make any more changes in this section? (y/n) 


n 


(eg. 


temperature [K] 
300.00 
400.00 
600.00 
700.00 
750.00 


0.,300. 


RETURN EN) 
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To change the input in any of the following sections,enter the item number 


of your Jchnoice 


1)NO MORE CHANGES 


2) CHANGE 
3) CHANGE 
4) CHANGE 
5) CHANGE 
6) CHANGE 
7) CHANGE 
&) CHANGE 
9) CHANGE 


PROPERTIES OF GLASS 
GEOMETRY OF THE WINDOW 
COEFFICIENTS 

FLAME RADIATION 

GAS TEMPERATURE 

HEAT TRANSFER COEFF. ON 
EMISSIVITY SOF SHOT LAYER 
NUMERICAL PARAMETERS 


HOT LAYER SIDE 


Would you like to save your input data in a file?(y/n) 


Would you like to save the output in a file? (y/n) 
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5. OUTPUT 


The output appears either on the screen in the form of numbers and can be written to 
files if desired by the user. File names must be provided by the user in order to save either 
Output. The file name must be in DOS format. The program deletes any pre-existing files 


with the same name. Thus to save all the runs, different file names must be used for each run. 


When saved in a file, the output begins with a listing of all the input vanables. Proper- 
ties of glass appear first, followed by geometry, coefficients, flame radiation, gas temperature, 
heat transfer coefficient on hot layer side, emissivity of hot layer and numerical parameters. 
Following the input variables is a table of the ume, the temperature of the side exposed to fire, 


the temperature of the side away from the fire, the dimensionless average temperature obtained 


approximately as 0,,. = £ (0(0)+8(1))*erf (Vt), and dimensionless time. On the screen, there 


is no listing of the input variables and the output begins with the table shown below in the 
example. After the table, the breakage time or a statement that breakage did not occur within 
the maximum mun time is given. Below this statement, the definitions of dimensionless time 
and dimensionless temperature and the values of the characteristic time, characteristic tempera- 
ture and the compression factor, g appear along with the values of average initial temperature, 
the average temperature increase at breakage and the average temperature at breakage. The 


example output that follows is for the default run. 


PHYSICAL AND MECHANICAL PROPERTIES OF GLASS 


i.Thermal conductivity [W/mK]= 0.7600e+00 

2.Thermal diffusivity [m*2/s]= 0.3600e-06 

3.Absorption length [m]= 0.1000e-02 

4.Breaking stress [N/m*2]= 0.4700e+08 

5.Youngs modulus [(N/m*2]= 0.7000e+11 

6.Linear coefficient of expansion [/deg C]= 0.9500e-05 


GEOMETRY 


i.Glass thickness [m]= 0.0063 
2.Shading thickness (mJ= 0.0150. 
3.Half-width [m]= 0.5000 


COEFFICIENTS 

1.Heat transfer coeff, unexposed [W/m*2-K]= 10.00 
2.Ambient temp, unexposed [K]= 300.0 

S., Emissivity *om glass = 99), 00 

4.Emissivity of ambient (unexposed) = 1.00 

FLAME RADIATION 


Number of points used for flux input: 2 


point # time [s] flux [W/m*2] 
aL OF OF 
2 100000 OR 


GAS TEMPERATURE 


Number of points used for temperature input: 17 


point # time [s] temperature [K] 
af O% 300.00 
Zz 2 ORO 303330 
3 20.00 503-96 
4 S000 304.60 
S ~ 40.00 305,57 
6 50.00 306.65 
7 60.00 308.20 
8 70.00 310.09 
9 80.00 342.64 

10 90.00 Sio.0.0 
sig 100.00 3197.56 
By: 110.00 325.19 
3 20. 00 3 Lees 
14 140.00 S53 271 
15 160.00 ae oS 
16 180.00 437597 
17 200.00 831.07 


HEAT TRANSFER COEFF. ON HOT LAYER SIDE 


Number of points used for heat transfer coeff input: 


point # time [s] h2 [W/m*2-K] 
" Os 50.00 
2 1000.00 50.00 


EMISSIVITY OF HOT LAYER 


Number, of ypoints! used for emissivity input: 2 
point # time [s] emissivity 

a QO. z.00 
74 1000.00 4,200 


NUMERICAL PARAMETERS 


1.Mazximum fractional error in soln= 0.000100 
2.5220 @ftime > ster "{'st="" 17000 

3.Mazimum run time [s]= 250.00 

4.Time interval for output [s]= 10.00 
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TEMPERATURE HISTORY 


Time Exposed Unezposed Theta Tau 

(s) T (K) THK) (Average) 

0. 300.0 300.0 QO. QO. 
10.0 30023 300.0 0.001 0.089 
20.0 300.6 300.0 0.002 04.279 
30.0 300.8 300.1 0.004 0.268 
40.0 SO Lak 300.2 0.006 OF 3I57 
50:70 301.4 300.4 0.008 0.446 
60.0 eOlne 300.5 0.012 07536 
70.0 802% B00? 0.016 ONG25 
80.0 302.9 B00 20 0.022 Orsi a4 
90.0 303... 7 BOdlsre 0.028 0.804 

100..0 304.6 301°. 6 0.036 02 8:9'3 
1.0)..0 305129 30 22.10 0.047 0.982 
120.0 30775 SU 0 0.061 07k 
230 0 309 27, 303.3 0.080 Le On 
140.0 feta 4° 304.2 0. 1S toe U 
15050 “316.4 305.4 0. 138 Goo s 
160.0 3.2 30740 0.180 1.428 
110.0 Sa se 3097.0 0.235 ih ote Ep 
180.0 S340 S3LLa6 0.304 L607 
L900 366.5 314.6 OwS36 e696 
200.0 434.9 S20¢S 1.033 on ao 
201.0 442.7 Bet. L082 12095 

Window breaks at time = 201 OOo ats) 

tau = t/tc, te = 112.0/'s,eAvg. theta = (Tav-Ta)7%ec, To = 7087 & 


g = 1.039 
Avg. T init = 300.0 K;¥Avg. Delta T = 77.2)K, Avg. T break?= Si7aeee 


i) 
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